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Abstract 



A theoretical framework is developed to describe experiments on the 
structure of epitaxial thin films, particularly niobium on sapphire. We ex- 
tend the hypothesis of dynamical scaling to apply to the structure of thin 
films from its conventional application to simple surfaces. We then present 
a phenomenological continuum theory that provides a good description of 
the observed scattering and the measured exponents. Finally the results of 
experiment and theory are compared. 
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Although there has been great interest in the theoretical modelling of 
molecular beam epitaxy (MBE) growth [1], there have been relatively few 
attempts to compare the detailed predictions of the models with experimen- 
tal results. Also, models developed from statistical mechanics have mostly 
concentrated on the shape of the growing interface rather than on the struc- 
ture of the grown material, whereas many experiments provide information 
about the whole film. Microscopic models of thin films have been developed 
[2] and successfully describe the structures of very thin films consisting of 
only a few atomic layers, but far less has been achieved in understanding the 
detailed structure of thicker films, 100-10000 A, that are of importance for 
many applications. In this paper, we present a dynamical scaling hypothesis 
for the structure of thin films involving the displacement of the atoms from 
their ideal positions. Based on this, we compare the results of x-ray scatter- 
ing studies of thin films [3,4] with the predictions of a continuum model for 
MBE growth developed by Villain [5]. 

It is well known [6], that for x-ray scattering from a single interface, 
the relative intensity of the specular reflection yields information about the 
mean-square width of the interface, < /i(x, t)^ >, in the direction perpendic- 
ular to the interface, while the form of the non-specular scattering provides 
information about the correlation function < /i(x, t)/i(x + r, t) >, where x 
and r are two-dimensional position vectors in the plane of the interface and 
the averages are taken over all positions x on the interface. The scattering 
from thin films requires more discussion [7]. For a wave vector transfer Q 
the scattering from a monatomic film is given by: 

/(Q) oc 5Zexp(i[Q||.(rj - rk) - Zk)]) (1) 

where (r^, zj) is the three-dimensional position of the atom j, the summations 
are over all the atoms in the film and the growth direction is parallel to 
the z-axis. We shall consider the scattering near the Bragg refiection, with 
Q = G parallel to the growth direction of the film so that Q\\ — q and 
Qz — G + Qz, with both q and being small compared with the reciprocal 
lattice vector G. The positions of the atoms in the film are written in terms 
of the atomic positions of the crystallographically perfect structure (Rj, Zj); 
the displacements from these positions for the real structure are written as 
Zj = Zj + Uj. Eq. (1) can then be rewritten as: 
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^(Q) = Eexp(^[q•(Rj - Rk) + qz{Zj - Zk) + G{uj - Uk)]) > (2) 



where it has been assumed that q and the displacements Uj are sufficiently 
small that terms of second order in these variables can be neglected. The sum 
over the z-components is from the fiat interface with the substrate Z — 0, 
to the height of the film Zj = d + hj, where d is the average thickness of 
the film and the hj are the fluctuations in the thickness whose average is 
zero. The Bragg-like component of the scattering Ib occurs for |q| = 0; 
if the distribution of the displacements and of the thickness fluctuations is 
Gaussian, Ib is given by: 

Ib{Q) = D{q,)5{q) exp(-G'2 <u^>), (3) 

where D{qz) is the Debye- Waller- like term arising from the fluctuations in 
the thickness of the film, given in turn by: 

D{q^) = (1 -2cos{q,d) exp{-ql < h'^ > /2) + exp{-ql < h'^ >))/ql 
This equation can then be further approximated to give: 

/^(Q) = 4g;2 sin2(g,rf/2)5(q) cxp(-G'' < > -ql < >) (4) 

The leading term in the diffuse or non-Bragg-like scattering is obtained 
by expanding the exponential term in Eq. 2 and assuming that the dis- 
placements from the ideal positions and the fluctuations in the heights are 
sufficiently small that only the leading terms need to be considered. The 
expression for g^^ = is then: 

Id{GI) = J1< {d^G^uJu]^ + hjhk) exp(iq.(Rj - Rk)) > (5) 

The displacements are averaged through the thickness of the fllm and are 
written u], so that the summations over j and k are now over all the lattice 
positions in the plane of the fllm. 

Equations 3 and 4 generalise to a thin fllm the expressions for the scat- 
tering from a single interface. The thin fllm has a conformal structure if the 
fluctuations of u are proportional to those of h. This corresponds to the case 
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when large fluctuations in the displacements are correlated with large fluc- 
tuations in the film thickness. The two terms in the scattering can then be 
combined into a single term. This assumption of conformality enables us to 
propose an extended dynamical scaling hypothesis for thin films. The scatter- 
ing is then solely a function of the fluctuations in the thickness/height instead 
of the two variables of the height and the atomic displacements. Equations 
3 and 4 form the basis of the scaling hypothesis for the scattering from thin 
fllms. 

With this assumption we now use conventional scaling relations for the 
thickness of the fllm (which we have argued is conformal with the atomic 
displacements) and their relationship to the scattered intensity, to deflne the 
critical exponents a and (3: 

u^{t) =< /i(x, tf >^ t^^, t <^L\t,L ^ oo (6) 
where L is the size of the fllm and z = a/ j3is the dynamical critical exponent, 

< /i(x,t)/i(x-Fr,t) >?a |r|^", |r| ^ oo (7) 

where the equalities above are valid asymptotically, and the scattered inten- 
sity is the Fourier transform of Eq. (7): 

/(q) Mq,^)M-q,^) > 

A continuum model of a surface was proposed by Villain [5] based on the 
Edwards- Wilkinson (EW) equation [8] with the addition of a surface diffusion 
term. The model leads to the differential equation: 

(9/i(x, t)/dt = -BV^h{^, t) + AV^/i(x, t) + ?7(x, t) (8) 

The flrst term on the RHS models the surface diffusion, the second is 
the normal diffusive term and ?7(x,t) is the noise that is assumed as usual to 
have a Gaussian distribution. When the surface diffusion term is neglected, 
the equation reduces to the EW equation, which has been used as an ap- 
proximation to the KPZ [9] equation to model the growth of multilayers [10]. 
We are unaware however of any comparison of the predictions of Eq. 8 with 
experimental data and in particular of the justiflcation of the presence of the 
different terms therein in the context of a speciflc experiment. 
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In his important work on the theory of MBE growth, Villain [5] pointed 
out the physical relevance of the terms in Eq. (8). In the presence of desorp- 
tion, the rate of change of /i(x, t) can be written as the difference between the 
chemical potential /i. at a point (x, t) on the surface and the average chemical 
potential /ic of the atoms in the molecular beam: arguing further that the 
consequent growth cannot depend on the orientation of the plane leads to: 

oc (/.(x, t) - /.e) = BV'h{^, t) (9) 

A further source of the V^h(x,t) term has been argued to be the downward 
funnelling of atoms to the lowest 'valley' available [11] to them. 

Turning to the surface diffusion term V^h(x,t), the rate of change of h is 
argued [4] to obey a continuity equation with respect to a current density j, 
which is itself the gradient of the chemical potential V^h(x,t): 

— ^ = -V.j(x, t) = -AVV(x, t) = -AV'hi^, t) (10) 

This term is characteristic of MBE growth, and is typically known as the 
noiseless Wolf- Villain equation [12,13]. It is a crucial ingredient in any serious 
model of MBE growth. Finally we remark that Eq. (8) is applicable to the 
motion of atoms on stepped surfaces [5]: additionally, downward funnelling, 
finite desorption and surface diffusion are valid physical processes and hence 
Eq. (8) is a reasonable model for describing the experimental conditions for 
MBE growth. 

Equation (8) is hnear and can be solved [5] exactly. If the system is also 
isotropic, the solution depends on g = |q|: 



{hiq, t)hi-q, t)) = {rjix, t)r)ix', t')){l - exp[-2( V + Bq^)t]}/{Aq' + Bq'), 
{h{q, t)h{-q, t)) = Co[l - exp(-27(g)t)]/7(g) 

(11) 

where Cq is a constant depending on the amplitude of the noise term and 
7(g) = Aq"^ -\- Bq^. If we assume that the film thickness d is proportional to 
time, as is assumed in theories of kinetic roughening [1] , the diffuse scattering 
is given by: 

lD{q) = C,[l - expi-2C2l{q)d)]/-f{q) (12) 
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where this expression follows from our dynamical scaling hypothesis for thin 
films, and Ci and C2 are constants. 

Figure 1 shows the x-ray scattering observed [3] from one of the thin films 
of niobium grown on sapphire single crystal substrates. The thickness of all 
the films is larger than the critical thickness and so the lattice parameters of 
the films are relaxing towards the bulk niobium lattice parameter. The solid 
line shows the result of a least-squares fit of the parameters Ci, C2, A and B 
in Eq. (12) to the data and we see immediately that Eq. (12) provides a very 
satisfactory description of the diffuse scattering. In particular for large q, the 
intensity is inversely proportional to 7(g) and so describes the wings of 
the scattering [3] very well, while at small q the theory has a finite value of 
the diffuse scattering as found in the experiments. Similar good agreement 
was obtained for the other layers and for the scattering around the (110) 
Bragg refiection. 

The full width at half maximum (FWHM) of the diffuse scattering in q 
was found [3] to vary with the layer thickness d as c^-o-5i±o o5 Assuming as 
usual that thickness is proportional to time, this suggests that the FWHM 
varies as t^^^"^. The correlation length of the model is defined [1] by ^(t) ~ t^/^ 
and the FWHM is proportional to 1/.^. The value of the exponent then 
imphes that the behaviour of the FWHM is dominated by the diffusive, 
V^/i(x,f), term in eq. (8). Since the line-shape for large q requires the 
quartic term, these two results strongly suggest that both the quadratic and 
quartic terms in Eq. (8) are needed to explain the experimental results. 

The experiments also measured the mean square fiuctuation of the struc- 
ture as a function of the thickness and found that u'^{d) ~ (^o-68±o.o5^ implying 
that the temporal roughening exponent is P — 0.34 ± 0.03.This value of /3 
is the value that is obtained ii A ~ for the solution of Eq. (8) for a one- 
dimensional surface; it is the value corresponding to the Wolf- Villain model 
i.e., with only the quartic term in Eq. (8). We therefore note that the be- 
haviour of both u{d) and the FWHM are in agreement with respectively the 
A — (quartic term only) and B — solutions (quadratic term only) of Eq. 
(8) for a one-dimensional surface. Normally we would expect that a two- 
dimensional system would exhibit logarithmic roughening [1], and we shall 
discuss below possible reasons for a dimensionality shift. 

We note that the spatial roughening exponent for the one-dimensional 
Wolf- Villain equation, (10), is a = 1.5 which is an indication of anomalous 
scaling behaviour [1]. Since the f3 exponent deduced from the measurements 
of u{d) is in good accord with the Wolf- Villain model we consider that at 
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least in some regime the films should exhibit anomalous scaling. There are 
some indications of this in Fig. 4 of ref. [3]; the scattering widths as a func- 
tion of i.e. perpendicular to the growth direction showed an interesting 
behaviour. For small qx (and corresponding to a lengthscale larger than the 
film thickness d), they were largely independent of qx, but for larger q^i they 
increased approximately linearly with q^. This was earlier [3] given the inter- 
pretation that in the latter case, the displacements with these wavevectors 
did not propagate through the film, while in the former case, the scattering 
profile was characteristic of the whole thickness of the film. Another way of 
seeing this, however, is to say that the distribution of the atomic displace- 
ments < ul > depends on the length scale over which it is examined. If 
the film is conformal, it is conceivable that this inhomogeneity in the distri- 
bution of atomic displacements should be reflected in surface fluctuations, 
which would then exhibit the same lengthscale dependence, leading in turn 
to anomalous scaling. 

Other experimental measurements of the spatial roughening exponent in 
metallic materials have given values of a = 0.79 ± 0.05 [14], 0.65 ± 0.03 [15] 
and 0.85 ±0.05 [16]. These values are intermediate between the values of the 
exponents predicted by the one- dimensional EW equation (V^/i), 1.5, and 
that of the Wolf- Villain model (V^/i), 0.5, which are the two limiting cases 
of our Eq.(8). Further investigations are needed to determine whether the 
experimental systems need entirely different theoretical models to describe 
them, or whether they could be treated as intermediate cases that could be 
accommodated within the general framework of our present work. 

Apart from the values of the exponents, the fits to the experimental re- 
sults provide information about the values of the parameters in Eq. (11). 
The measurements were made in such a way that the scattered intensity 
could not be reliably scaled from one sample to another and so the constants 
A and B were obtained by adjusting the overall scale factor to give the ap- 
propriate total scattered intensity for the thickness of the film. The results 
for parameter A are shown in Fig. 2. There is considerable scatter in the 
results at least partly because the parameters in the least squares fits are 
strongly coupled, and both A and B are consistently larger when deduced 
from the (110) data in comparison with the (220) data. As discussed in [3], 
some of the differences between the films may arise because they were grown 
at different times and capped in different ways. The difference between the 
results obtained from scattering near the (110) and (220) Bragg refiections 
may arise from a failure of the small displacement approximation. Also, the 
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small value of A obtained for the film with a thickness of 1500A may be be- 
cause the correlation length is then so small that it is difficult to distinguish 
between the Bragg component and the diffuse scattering. 

Figure 3 shows the thickness dependence of the B parameter. This term 
arises from the quartic term in the differential equation (8) and so to compare 
it with the A term we need to scale it by the square of a length. In Fig. 3 
we show B/(f because we found empirically that this quantity was largely 
independent of d. 

A full microscopic theory of the structure of thin films is not yet available. 
The distortions from a perfect structure arise presumably from the existence 
of threading dislocations as explained, for example, by Kaganer et al [17]. 
Although their theory can be used to describe thin films with very low dis- 
location densities or for very thick films, it is the intermediate region that is 
important for many applications. 

In conclusion, we have proposed a scaling ansatz to describe the devel- 
opment of the structural deformations in thin films and have applied it to 
describe the structure of Nb thin films grown on sapphire substrates. We 
note that other films, both semiconductors and metals give similar experi- 
mental results. The scattering is well described by the profile proposed by 
Eqs. (8, 12). The exponent measured for the spatial correlation length shows 
that it is dominated by the diffusive interactions (quadratic term) in the phe- 
nomenological differential equation (8), while the quartic term in Eq. (8) is 
needed to explain the temporal roughening exponent. We believe we may 
in this sense have provided a minimal model to explain various quantitative 
features of the experiment in ref . [3] . 

Normally we would expect that a two-dimensional system such as the 
corresponding to the experiment in [3] would exhibit logarithmic roughening 
[1], rather than the power-law roughening indicated by our measurement 
of the temporal roughening exponent, (3. Since the simplest theories state 
that two-dimensional surfaces are at the critical dimension for which the 
roughening shows logarithmic behaviour when the diffusive term is present, 
the power law result is unexpected. We note that good agreement with the 
observed exponent is obtained instead for the linear one-dimensional theory 
presented here. 

The most likely explanation for this is that since the measurements con- 
cerned were the result of synchrotron experiments, there was a high degree 
of anisotropy in their resolution. In particular, in the experiments [3], the 
experimental resolution element perpendicular to the scattering plane was 
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much larger than in the other directions; the integrated data over the x-y 
plane could also have contained some highly anisotropic averaging. It could 
also be that there is in the experiment [3] a dimensionality shift analogous to 
that found in the random field model [18] due to the initial defects present 
in the substrate. Another remark concerns the fact that the ratio B/dP is 
approximately constant (Fig. 3), while A decreases with increasing d, lead- 
ing to the increasing dominance of B for large film thicknesses. This could 
be because as the thickness of the film increases, with its length always re- 
maining constant, the approach to asymptoticity in eq. (8) is slowed down, 
as a result of which the higher-order nonlinear term (whose coefficient is S, 
in this case) appears to win out. We hope that the relative success of our 
minimal model in explaining quantitative observations on the sapphire-Nb 
film analysed in [3] will lead to some of these issues being resolved, and to 
new theoretical developments in the study of thin films and multilayers , for 
which the thickness is in the important range between 100 and 10000 A. 
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Figure captions: 
Fig.l: 

X-ray scattering intensity from a transverse scan of the wave- vector transfer 
through the (220) Bragg refiection from the Nb layer, 20nm thick. The 
observed intensity is fitted to a Gaussian to represent the Bragg reflection 
and the function proposed in the text for the diffuse scattering. 
Fig. 2: 

The parameter A as a function of the layer thickness as deduced from scat- 
tering near the (110) Bragg reflection (flUed points) and the (220) reflection 
(open points). 
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Fig. 3: 

The parameter B/(P as a function of layer thickness. The data were obtained 
from scattering near the (110) Bragg reflection (filled points) and the (220) 
reflection (open points). 



References: 

1. J. Krug, Advances in Physics 46 139 (1997) 

2. S. Das Sarma and P. Tamborenea Phys. Rev. Lett. 66 325 (1991) ;S. 
Das Sarma and S. V. Ghaisas, Phys. Rev. Lett. 69 3762 (1992) 

3. A. R. Wildes et al., J. Phys: Condens. Matt. 10 L631 (1998) 

4. P. M. Reiner et al., Phys. Rev. 45 11426 (1992), G. Gutekunst et al., 
Phil Mag. A 75 1329, 1357 (1997) and P. F. MiceUi et al, Physica 221 230 
(1996) 

5. J. Villain, J. de Physique I 1 19 (1991) 

6. S. K. Sinha et al., Phys. Rev. B 38 2297 (1988) 

7. V. Holy, U. Pietsch, and T. Baumbach; High resolution x-ray 
scattering from crystalline thin films, Springer Tracts in Modern Physics, 

Vol.149, Springer- Verlag (1999). 

8. S. F. Edwards and D. R. Wilkinson, Proc. Roy. Soc. A 381 17 (1982) 

9. M. Kardar et al., Phys. Rev. Lett. 56 889 (1986) 

10. T. Salditt et al., Phys. Rev. Lett. 73 2228 (1994) 

11. D. E. Wolf and J. Villain, Europhys. Lett. 13 389 (1990) 

12. J. W. Evans, Phys. Rev. B 43 3897 (1991) 

13. G. L. Zhou and C. P. Flynn, J. Phys: Condens. Matt. 9 L671 (1997) 

14. Y. L. He et al., Phys. Rev. Lett. 69 3770 (1992) 

15. X. Yan and T. Egami, Phys. Rev. B 47 2362 (1993) 

16. P. P. Swaddhng et al., Phys. Rev. Lett. 73 2232 (1994) 

17. V. M. Kaganer et al., Phys. Rev. B 55 1793 (1997) 

18. Y. Imry and S. K. Ma Phys. Rev. Lett. 35 1399 (1975) 



10 



